Infinite families of superintegrable systems 
separable in subgroup coordinates 



Daniel Levesque, Sarah Post and Pavel Winternitz 

Abstract. A method is presented that makes it possible to embed a subgroup 
separable superintegrable system into an infinite family of systems that are integrable 
and exactly-solvable. It is shown that in two dimensional Euclidean or pseudo- 
Euclidean spaces the method also preserves superintegrability. Two infinite families 
of classical and quantum superintegrable systems are obtained in two-dimensional 
pseudo-Euclidean space whose classical trajectories and quantum eigenfunctions are 
investigated. In particular, the wave-functions are expressed in terms of Laguerre and 
generalized Bessel polynomials. 



1. Introduction 

The purpose of this paper is to show how infinite families of superintegrable systems 
are obtained from individual known ones which allow separation of variables in at least 
one subgroup type coordinate system. In this paper we restrict to two-dimensional 
spaces that are real forms of the complex Euclidean space ^(C). The systems obtained 
in real Euclidean space are already well known [5j [161 HH H31 H5] . We also obtain two 
infinite families of superintegrable and exactly-solvable systems on the pseudo-Euclidean 
plane, E^i, previously obtained in [17]. These systems are related by coupling constant 
metamorphosis and we investigate their classical and quantum systems in depth. 

We recall that a superintegrable system is one that allows more integrals of motion 
than degrees of freedom. The number of functionally independent classical integrals in 
an n-dimensional space is at most 2n — 1. A superintegrable system in two-dimensions 
will allow precisely 3 independent integrals. 

Superintegrable systems are rare. As a matter of fact, Bertrand's Theorem [2j [8] 
tells us that the only spherically symmetric superintegrable potentials are V(r) = a/r 
or V(r) = ar 2 (in any dimension). On the other hand, these are the two most important 
potentials in physics. 

Superintegrable systems in classical mechanics are of interest mainly because, in 
these systems, all bounded trajectories are closed and the motion is periodic |36j. In 
quantum mechanics, the energy levels are degenerate and superintegrable systems have 
been conjectured to be exactly solvable [44J. We recall that in quantum mechanics 
a system is called exactly solvable if its entire energy spectrum can be calculated 
algebraically and its wave function are polynomials in the appropriate variables, 
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multiplied by an overall gauge factor (the ground state wave function) (9j IHJ 1^8] - 
Both in classical and quantum mechanics, the integrals of motion for superintegrable 
systems form interesting non-Abelian algebras (under Poisson or Lie commutations, 
respectively). The algebras are polynomial ones [IJ [TTJ [3Q1 132] that exceptionally simplify 
to Lie or Kac-Moody algebras [3]. 

Until recently, most studies of superintegrable systems concentrated on quadratic 
superintegrability, i.e. the case where all integrals of motion are at most second- 
order in the momenta, see e.g [7J [19J, EDJ EU 151] . Second-order superintegrability is 
related to multi-separability in the Hamilton- Jacobi and Schrodinger equations. Indeed, 
in real Euclidean spaces, superintegrability and multi-separability are in one-to-one 
correspondence. 

It was shown in [7] that exactly four families of quadratically superintegrable 
systems exist in E 2 (M.). One of them allows separation of variables in two subgroup 
type coordinates, namely 



Hi = p 2 + -p 2 e + u 2 r 2 + - — + „ (1.1) 



» 1 2 i 2 2 i a @ 

" ~^Pe + ^ r + — — + — — 

2 i 2 i 2/2, 2\ , a , @ 
Px+Py + U [X +V ) + -o + ~2 



(polar and Cartesian coordinates), two in one subgroup and one non-subgroup type, 
namely 



a 



H 2 =p 2 x +p 2 v + co 2 (x 2 + Ay 2 ) +Py + - (1.2) 



x 



(cartesian and parabolic) and 



TT 2 , 1 2 . a . 1 / P 1 \ f-l Q\ 

H 3 = P r + -^Pe + - + ~2 — z~a + ' L3 

r 2 r r 2 \cos 2 sin 9 J 

(polar and parabolic). The fourth system allows separation in two different parabolic 
systems and is not relevant to this article. 



2. Adding a parameter to integrable systems that separate in 
subgroup-type coordinates 

Let us first consider the two-dimensional Euclidean space E 2 (C) with its isometry group 
E(2, C) and its Lie algebra e(2, C) ~ {L 3 , Pi, P 2 }. This algebra has three non-equivalent 
one-dimensional subalgebras, {-Pi}, {Pi+iP 2 } and {L 3 }. Diagonalizing Pi, Pi+iP 2 or L 3 
simultaneously with the Laplace operator A = P 2 + P 2 corresponds to the introduction 
of subgroup coordinates, namely cartesian {u, v}, light-cone {£ = u + v , r\ = u — v} and 
polar {p = \Ju 2 + v 2 , 9 = arctan^} coordinates, respectively. 

A classical Hamiltonian on E 2 (C) separable in orthogonal subgroup type 
coordinates (i.e. Cartesian or polar coordinates) will have the following form 

H = p\ + fMi) + Hqi)X, X = p\ + f 2 (q 2 ), (2.1) 
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where (gi, q 2 ), (pi,P2) are the classical coordinates and the conjugate momenta on four- 
dimensional phase space. In these coordinates, the kinetic energy, or free Hamiltonian, 
is 

A=pl + ^( qi )pl (2.2) 

and V = /i(<7i) + ^(gO/^fe) is the potential. In E 2 (C), the only two inequivalent 
choices for ip are 

1 Cartesian coordinates 
1/qf polar coordinates 



(2.3) 



There is a natural quantized version of (2.1) obtained by taking (91,52) as the position 



operator and (pi,P2) as the momentum one, though there may be quantum corrections 



to (2.1). 



We observe immediately that system (2.1) is integrable and that X is an integral 
of motion. Further, if we introduce a parameter k by replacing X — y k 2 X, the obtained 
system remains integrable but the metric associated with the kinetic energy is changed. 
The system is mapped back into flat space by making an additional change of variables 
q 2 — y kq 2 , p 2 (1/&)P2- The new Hamiltonian and integral of motion are 

H = Pi + /i(<7i) + i!>(qi)X, X = p 2 2 + k 2 f 2 (kq 2 ). (2.4) 

The kinetic energy is the same as in ( [2.1 ), however, the potential term f 2 has acquired 
a parameter k. The Hamiltonian (2.4) represents an infinite family of integrable 



Hamiltonians which are separable in the original coordinate system (qi,q 2 ). 

Let us know apply the above procedure to known two-dimensional superintegrable 



systems. Starting from the system (1.1) in polar coordinates, we arrive at the TTW 
system 

a (3 



H 



TTW 



1 k 2 

2 1 2 1 2 2,^ 

Pr + - Pe + ur + — 



cos 2 k9 + sin 2 k6\ ' ^'^ 
It is integrable and exactly solvable and was conjectured to be superintegrable for 
rational k [1SI US]- Several authors then proved that (2.5) is superintegrable for all 



rational k, both in classical and quantum mechanics [Ml EHJ H2] (without any quantum 
corrections). 



The system (1.1) is also separable in Cartesian coordinates and the same procedure 



applied to both x and y leads to the anisotropic singular (or caged) harmonic oscillator 

urn 



Pt + P 2 y + " 2 (ktx 2 + k$y^ 



ft 



+ 



/3 



x- 



y 



2 ■ 



(2.6) 



It is clear that this procedure applied to (1.2) in cartesian coordinates will give a subcase 



of (|2.6|) and in fact (|1.2|) is a member of the infinite family (|2.6|). 

f 

f3k 2 



Finally, starting from the deformed Coulomb system (1.3), we arrive at the system 

ak 2 



Hi 



PW 



K 
+ — + 



r 2 cos 2 kO ~^ r 2 sin 2 k9 



(2.7) 
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which was proven |JT] to be superintegrable and equivalent via coupling constant 
metamorphosis to the TTW system (2.5) (with k — > 2k). Thus, (2.5), (2.6) and (2.7) 
are the only infinite families of superintegrable systems, obtained in this manner, in 
E 2 (C) that also exist in the real Euclidean space E 2 = E 2 (M). 

Two further superintegrable systems exist in E 2 (C) that are separable in one 
orthogonal subgroup type system of coordinates (polar coordinates) and one non- 
subgroup type (hyperbolic coordinates). The corresponding Hamiltonians in Cartesian 
coordinates are 

H^. A+A + ^ xl + A)+aj ^ i+P ^l^ (2 . 8) 

and 

H DC = p\+p 2 2 + Jj— + a \ + d } —.(2.9) 

^x\ + x 2 2 (xi + ix 2 y ^x\+xl{xi-ix 2 ) 

These two systems were identified as E7 and E17 in [25J. They separate in polar 
coordinates on E 2 (C) 

Xi + ix 2 = re* 9 , Xi — ix 2 = re~ w , 

and we have 

H DO =pl + ±- 2 pl + co 2 r 2 + 1 {ae ue + Pe 2W ) , (2.10) 

H DC = p 2 r + \p 2 + - + \ (ae™ + 0e i6 ) . (2.11) 

In real Euclidean space, H DO and H DC would be respectively a deformed harmonic 
oscillator and a deformed Coulomb potential and the deforming potential in both cases 
is complex. 

Still in complex Euclidean space, let us introduce a parameter putting 6 — > k9. This 
preserves one second order integral and separation of variables in polar coordinates. 
Following the procedure outlined above with a scaling of the coupling constants, we 
obtain the integrable families [?] 

H?° =p 2 r + ±p 2 e + co 2 r 2 + I (ae Alkd + /3e 2 ^) (2.12) 

H™ =pl + ^P 2 e + * + ^ (ae™ + 0e™) . (2.13) 
The two Hamiltonians H®° and are related by coupling constant metamorphosis 



T5l [271 129] . In fact, if we consider the Hamilton- Jacobi for if 



DO 

k 



Pr + ^Pe + "V + y 2 (^ 4lk9 + fc 2 ™) = E, (2.14) 



and solve the equation for u, we obtain 



1 ( PI + y 2 pi) - % + ^ + ^) = "- 2 - (2-15) 



.2 
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The change of variables p = r 2 /2, 



29, changes the left-hand side of (2.15) 



into the Hamiltonian (2.13) in the new polar coordinates p, with K = —E/2 and the 



constants a, (3 scaled by 1/4. 

Just as in the previous known examples of infinite families of Hamiltonians, there is 
no guarantee that the infinite family is superintegrable. The sup erintegr ability of each 
of the families must be verified. In the following sections, we prove that the system 



(2.12) is superintegrable in both the classical and quantum case and hence so is (2.13) 



as a result of the coupling constant metamorphosis. 



3. Infinite families of superintegrable systems in the pseudo-Euclidean 
plane 



Let us now restrict to the real pseudo-Euclidean plane Ei t i and study the system (2.12) 
In cartesian coordinates, we have 



H?° =pl+p 2 2 + co 2 {xl + x 2 2 ) + 



rp 2 I rp 2 



a 



xi + ix 2 
x\ — 1x2 



2k 



+ /3 



X\ + IX2 
X\ — IX2 



Let us put 

X\ = U, 1X2 = V, 11,1)6 

We obtain a Hamiltonian in En namely 



(3.1) 



(3.2) 



tt 2 2,2/2 2\ 1 

Hk =p u -p v + uj [u -v ) + 



1 



a 



u- 



U + V 



u — V 



2k 



+ P 



U + V 



u — V 




a,pe R.(3.3) 



Without loss of generality, we can assume that k > and either a = /J = 0ora^0; 
indeed a = is equivalent to a 7^ 0, /3 = with k — > k/2. In the former case, u, v are 
allowed to take all real values. In the latter case, we require that u 2 — v 2 be positive 
unless k is a positive integer. If A; is a positive integer, the requirements are relaxed so 
that it is only required that u — v 7^ 0. 



The Hamiltonian (3.3) allows separation of variables in pseudo-polar coordinates 
r cosh 7, v = r sinh 7, 



u 



r 2 = u 2 —v 2 , tanli7 



v 
u 



with r G K for the "time-like" region of the pseudo- Euclidean plane. In the "space- 
like" region, we must interchange u and v. An equivalent but more convenient set of 
coordinates in are modified pseudo-polar coordinates 

u + v a + 1 fp a — 1 jp 

V = 2 V a' 



P 



2 2 
u — v , 



a 



u 



u — V 

-00 < p < 00, —00 < cr < 00, 



2 
P 
a 



(3.4) 



(u-vY > 0. 



(3.5) 



These coordinates cover Ei t i with the line u — v = removed. The Hamiltonian in these 
coordinates is 

H k = App 2 p - A-pl + u 2 p+- (aa 2k + (5a k ) . (3.6) 
P P 
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We shall show below that the corresponding Hamiltonian system is superintegrable 
both in classical and quantum mechanics for all rational values of the constant k. The 



potential in (3.6) is real and finite for < p < oo, < a < oo and we restrict to this 



region when either a or f3 is different from zero. 



4. The classical system 

4-1. Separation of the Hamilton- Jacobi equation. Trajectories 



In this section, we solve for bounded trajectories of the Hamiltonian system (3.6) in the 
case that a^O. The case of a = (3 = is the harmonic oscillator in E^i, similar to the 
one in E 2 although the energy can take negative values. The trajectories can be directly 



obtained from separation of variables in the coordinates u, v as discussed in Appendix 

m 



Consider a bounded trajectory in the following regime, 
< Pi < p < P2, < <Ji < a < a 2 . 



We write the Hamilton- Jacobi equation for Hk as in (3.6) as 

2 



2k + po- h ) 



E, 



and separate variables, putting 

2 



4P 2 



dp 



i 2 2 
+ UJ p 



Ep = 4a 2 



d_ 



aa 



2k 



f3a k 



A 



with 



A e 



S = S 1 (p) + S 2 (a) - Et, 
The functions Si and S 2 are 

1 r p2 

Si(p) = -o / V-w 2 + Et- 1 + At-Ht + a, 
S 2 ( a ) = -- VAt- 2 + aT 2k - 2 + PT k - 2 dT + c 2 , 



p > 0. 
a > 0. 



To calculate the trajectories, we solve the equations 

dS _ dSi 

dE~ ~dE 
OS _ OS! dS 2 

dA~ ~dA + ~dA 



t = 5 u 



do 



(4.1) 



(4.2) 



(4.3) 



(4.4) 
(4.5) 

(4.6) 
(4.7) 



where Si and 5 2 are arbitrary real constants. Beginning with (4.6), we obtain 

1 f p2 dT , . 

t + 5 1 = -- . 4.8 

AJ p V-uj 2 t 2 + Et + A V ; 

Immediately from this equation, we see that unless u 2 > 0, the trajectories will be 
unbounded at t — > oo. Further, since the polynomial —u 2 p 2 + Ep + A is negative for 
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large r, the integrand (4.8) will take real values if and only if the polynomial has two 

(4.9) 



distinct real roots. With this requirement, the integrated form of (4.8) is 

-2co 2 p 2 + E' 



P 



2co 2 



E — \/ D\ sin ( — Au(t + 8i) + arcsin 



, D 1 



E 2 + 4u 2 A > 0. 



From these equations, we see that the limits on the values of p are 



E 



Pi 



< < 



P S P2 



2tu 2 ~ r ~ r 2tu 2 
and hence the equation for p simplifies as 
1 



P 



2tu 2 



E+ V^Di cos (4u(t + 5i)) 



(4.10) 



(4.11) 



Note that that p has a period of ir /2u in t. In the case that a = (3 = 0, we have the 
trajectories for the harmonic oscillator, compare with (A. 5). In this case A is positive 
and so pi < < pi- On the other hand, if we assume either a or is non-zero, this leads 
to the requirement p > which implies both E > and A < 0, as well as E 2 +4uj 2 A > 0. 

Next, we solve the equation (4.7) in the sector where p, o > 0. The first relevant 
integral is, 



with 



dS 1 
~dA 

dSx 



4 



P2 



r-V-w 2 + Et- 1 + At~ 2 



1 



(2Ap- 1 + E 

arcsin 



OA Ay/-A f~ ~~ V VDi 
The second relevant integral in (4.7) is, for k ^ 
dS 2 -1 n dT 
~dA ~ ~4 

1 



A < 0, 



(4.12) 



(4.13) 



T 



■^JAt- 2 + aT 2k ~ 2 + @t 



k-2 



4k 



<T 2 



—k\ 



\J Ai 2k + (3T~ k + a 



which is integrated as 



dS2 
dA 



Ak^A 
D 2 



arcsin 



(2Ao- k + (5 



/D 2 
2 AaA > 0, 



arcsin 



2Aa 2 k + f3 



ID* 



(4.14) 



(4.15) 



a 2 



^Dn 



2\A\ 



Thus, all bounded trajectories satisfy E > and A < with 
i — S o~i < a < a 2 < 



2\A 



2 A 



(3>0, a<0. 



(4.16) 



The condition D 2 > provides a fundamental restriction on the values of the separation 
constant, namely 

/3 2 



< -A< 



4a 



(4.17) 
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We are now in a position to calculate the trajectories. Let us define new variables 
_ 2Ap- 1 + E TTr 2Ao- k + (3 



W, 



and the constant 



T k = k (a5 2 V^A + |) + arcsin ( 



D 2 

2A{-(3 + ^D~ 2 ) + (3 



(4.18) 



(4.19) 



Putting (4.13) and (4.15) into (4.7), we obtain an equation for the trajectories, namely 

(4.20) 



arcsin W k + k arcsin Z = T k 



or equivalently 



a 



-k 



-1 

2A 



D 2 sin (—k arcsin Z + 



We already know from (4.11) that p and hence Z has period r 



2aJ 



(4.21) 
. Since arcsin is 



a multivalued function, equation (4.21) does not prove that a has the same period, or 
that it is periodic at all. On the other hand, if k is rational, equation (4.21) can be 
reduced to a single valued implicit function of p and a and this proves that the motion 
is periodic. Put 



k — p/q, p, q G N. 



(4.22) 



We rewrite (4.20) as 

cos C k = cos [q arccos Wk + p arccos Z] , 
Making use of the Chebyshev polynomials defined as, 



7T 



C k = (p + q)--qT k . (4.23) 



T n (x) = cos (n arccos(x)) , 



UJx] 



sin ((n + 1) arccos(a;)) 



sin arccos x 



(4.24) 



equation (4.23) becomes 



cos(C k ) = T q (W k )T p (Z) - U q ^(W k )U p ^(Z)^(Wl - 1)(Z 2 - 1). (4.25) 
Using the explicit formulas for the Chebyshev polynomials 



[n/2] 



^) = ]T(^- 2 V-i) J 



j=0 



[n/2] 

^) = E(^ n ~ 2 v-i) j 

j=0 



where n is a positive integer, we can now write (4.25) as 



[9/2] 



b/2] 



7=0 



7=0 



W-l)(^-l) 



[(g-i)/2] 



E 



(g-l)-2i m/2 
fc 



J'=0 



KP-l)/2] 

E (W^- 5 ^ 



Let us summarize the results on the trajectories. Without loss of generality, we 
have assumed k > and «^0. 
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(i) Bounded trajectories are obtained for 



> a,0 < P, 



< -A < 



E 2 



< -A < 



J 2 



(4.27) 



4w 2 ' ----- 2a 
Thus, for this system in contrary to the case of the real Euclidean plane 

E 2 [13 US], we obtain a fundamental upper bound on the value of the separation 



constant A (4.27) 



(ii) For rational k, all bounded trajectories are periodic. 
4-2. Superintegrability of the classical system for rational k 

The classical trajectories can be used to explicitly construct an integral of the motion 
which is polynomial in the momenta (see e.g. JT7J E3J SI])- The additional classical 
integral of the motion was first constructed in [17]. It is also possible to construct the 
integral directly from the bounded trajectories by replacing the conserved quantities A 
and E by their phase-space counterparts 

E -)• H, A ->■ Aa 2 p 2 a - aa 2k - 0a k . 

The relevant quantities are 



(4.28) 



Z = 

w k 



Ap 2 p 2 p + Aa 2 p 2 a + u 2 p 2 - ao 2k - f3a k 



g a 2-k p 2 



2aa k - /3 



Wo 



(1-Z*)(1-W 2 ) 



16a 1 k p p p a {aa 2k + (3a k - ^a 2 p 2 ) 



(4.29) 
(4.30) 
(4.31) 



D l = H 2 + 4qu 2 (4o- 2 pI - aa 2k - f3a k ) 2 (4.32) 
D 2 = (3 2 - Aa(Aa 2 p 2 a - aa 2k - (3a k ). (4.33) 

Using the parity properties of the Chebyshev polynomials, the following function on 
phase space is a constant of the motion which is polynomial in momenta of total degree 
no greater than 2p + 2q, 

g 



'Di V^2 



T g (W k )T p (Z) - U^iW^U^^JiW 2 - 1)(Z2 



(4.34) 



Here we see the explicit link between the integrals of motion and the bounded 
trajectories. As discussed elsewhere, see e.g. [26], this is not necessarily the integral 
of motion of lowest degree in the momenta. Nevertheless, this integral shows that the 



Hamiltonian (3.6) is classically superintegrable for rational k. 



4-3. Examples of trajectories 



We shall now calculate and plot the trajectories directly from (4.25). For k = p/q 

(4.35) 



integer, (4. 25) becomes 



cos(C p ) = W P T P (Z) - U P ^(Z)J(W^ - 1)(Z 2 - 1 
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k = 1 





\ 

) 


j -0.5 








jiT 0.5- 




(: 




k = 


= 2 




k = 3 



Figure 1. Trajectories of specific cases for k integer. The constants are chosen to be 
A = -1 , a = -2, (3 = 6, E = 20, u = 3 and <5 2 = tt/32. 




jfe = 1/3, 5 2 = 3tt/12 



jfe = 1/2, 5 2 = vr/6 



fe = 3/2, <5 2 = tt/12 



Figure 2. Trajectories of specific cases for rational k. The constants are chosen to be 
A = -3/2, a = -1,/3 = 3,E= 20 and oj = 4. 



and trajectories for k = 1, 2 and 3 satisfy: 

fc = 1 : cos^d) = WiZ - y/(W? - 1)(Z 2 - 1) 
k = 2 
k = 3 



cos(C 2 ) = W 2 (2Z 2 - 1) - 2Z V /(W 2 -1)(Z 2 -1) 



cos{C 3 ) = W 3 {4Z 3 - 3Z) - (4Z 3 - 1) ^/(Wf - 1){Z 2 - 1). 
The trajectories are plotted in figure 1 with constants chosen as indicated in accordance 



with (4.27) 



The trajectories for k = 1/3, 1/2 and 3/2 satisfy: 



k = § : cos(tf 1/8 ) = (4W* /S - 3Wi /a )Z - (4^% - lW(W? /3 - 1)(^ 2 - 1) 



k = \ : cos(C 1/2 ) = (2^ 2 /2 - l)Z - 2W 1/2 J(W? /2 - 1)(Z 2 - 1 



jfe 



3 . 



cos(C 3/2 ) = (2Wf /2 - 1)(4Z 3 - 3Z) - 2^ 3/2 (4Z 2 - l) A /(W| /2 - 1)(Z 2 - 1 



3/2 



The trajectories are plotted in figure 2 with constants chosen as indicated in accordance 
with ((4271) 
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In quantum mechanics, the Hamiltonian on pseudo-Euclidean space E(l, 1) is 
Hk 



d 2 d Aa 2 d 2 Aa d 2 1 , 2k Q k . 

-Ap— -4— + — — - H — — h up + - (aa 2k + (3a k ) 

op 1 op p oa 2 p aa p 



The integral associated with separation of variables is 



-Aa 



o 2 



d 



- Aa— - aa 2k - fia k . 
da 2 da 



Here again, we assume that k > and a is not zero. 
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(5.1) 



(5.2) 



5.1. Wave functions and energy spectrum 

We look for wavefunctions which are also eigenfunctions of the integral L\ and hence 
allow separation of variables 

H k fy = E~$, L 1 ^ = A^, V = R(p)S(a). (5.3) 

The functions R{p) and S(a) satisfy the following equations 

A 



' d 2 d 2 
Ap— +A—-u 2 p+- + E 
op z op p 



R(p) = 0, 



Aa' 



da 2 







+ Aa— + aa 2k + (5a k + A 
aa 



S(a) = 0, 



(5.4) 
(5.5) 



respectively. 



Equation (5.4) can be solved in terms of the confluent hypergeometric function and 



indeed we have two solutions 



1 + ^Ta 

2 



e 2 ^(a;p/2) ± ^ 1 F 1 ( t ± ^ - 1 ± V^A; 



up 



(5.6) 



that are linearly independent as long as 1 ± y— A is not an integer [10] . The requirement 
that the solution be continuous at the origin necessitates the choice of R + only. The 
requirement that R{p) be square integrable implies that the iiq must be a polynomial. 
This provides the quantization condition 

1 V^A _ E_ 

2 + 



-m, 



.e. E = 2u(2m + 1 + \f—A), m G Z. 



(5.7) 



2 4w 

The bound state eigenfunction R(p) can alternatively be expressed in terms of the 
Laguerre polynomials as 



R(p) 



e 2 



uip 



(up) ; - 



1 l: 



\up). 



(5. 



Equation (5.5) can also be reduced to that for the confluent hypergeometric series 
by putting 



A O 7„ 
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(the other solution is singular at a = 0). Square integrability again necessitates that 
the 1-F1 functions should reduce to a polynomial and we obtain a quantization condition 
for the eigenvalue A, namely \ — j-J=z + 

An = k 



2 4V = afc 
P 



2k 



-n i.e. 



2JU 



2n 



-a 



a < 0,/3 > 0. 



(5.10) 



Note that a and f3 satisfy the same bounds as in the classical case (4.27). We see that 
the quantized Hamiltonian (5.1) only allows a finite number of bound states. Indeed, 

(5.11) 



the condition y/—A n > implies 



n < 



P 



2 \2k, 



1 



-a 



In particular, the ground states exists only if 



P 



> 1. 



(5.12) 



2Ay-a 

In analogy with the classical case, the separation constant A can take values only in 
a finite range (in this case a quantized one). Let us define N as the largest integer 



satisfying (5.11) and define 



B 



P 



2kJ^i 



-2N-1. 



(5.13) 



We then have 

^A~ n = k[B + 2(N -n)}. (5.14) 
The functions S(cr) of ( |5.9 ) can be expressed in terms of the Laguerre polynomials as 

Sir,) = e-^ N -"L^ N -«\ V ), V = JE™!L. (5.15) 



k 



Changing variables in (5.5), we see that S{rj) satisfies a self-adjoint differential equation, 
namely 



— (v-s ) 

drj \ drj 



1 A 1 P 

—A + -7] + A1 

Akr] 4 Aky/^a 



S = 0. 



(5.16) 



It follows that the eigenfunctions S n corresponding to different values of n are mutually 
orthogonal 

S n S n ,-t = N n 5 n , n ,. 0<n,n'<N. (5.17) 
rj 

We mention that this orthogonality does not follow from the orthogonality of the 
Laguerre polynomials L° in (5.15) since the index a depends on the lower index n. 

On the other hand, the functions S(a) solving (5.5) are directly related to the 
generalized Bessel polynomials y n {x, a, b) introduced by Krall and Fink [28]. They satisfy 
the equation 

x 2 y" + (ax + b)y' = n[n + a — l)y, y — y n (x , a, b) . (5.18) 
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Indeed, let us parameterize a in (5.5) as 

P 2 



a 



4fc 2 (a-2) 
change variables putting 

k 1 



a ^2, 



n 



and express the solution of (|5.5|) as 

SJs) 



2e *>B n (s,a,l). 



(5.19) 



(5.20) 



(5.21; 



The polynomials B n (s,a,l) satisfy (5.18) and are hence the generalized Bessel 
polynomials defined as [281 EH] 



B n (s, a,b) = ) j C n> k (n + k 



2 )W 



x 



(5.22) 



fc=0 

where C„ h is the binomial coefficient and 



,(*) 



z(z — l)...(z — k + l) 



is a Pochhammer symbol. Note that the orthogonality of the Bessel polynomials over a 
finite range of indices is consistent with the results of [6]. 

5.2. Superintegrability of the quantum system for rational k 

Let us now use the recurrence relation approach of Kalnins, Kress and Miller [22] to 



show that the quantum Hamiltonian is superintegrable. The main results of section [5TT 
are that the energy and separation constant for the system are 



E„ 



2a; 2m - 2kn + 1 - k + 



k -2n + 



ft 



- 1 



2^ak 

and the wave functions are given by 
with 



S n ( V ) = d n e~ ri ^ 2 r]^ +N ~ n L^ +2N ~ 2n (i]), V 



—aa 



(5.23) 
(5.24) 

(5.25) 

(5.26) 
(5.27) 



For rational k — p/q, we can exploit the degeneracy of the energy with respect to the 
transformations 

m— >m + p, n— >n + q 
m —7- m — p, n — )■ n — q 

to show that such a system is superintegrable. The additional integral of motion will 
be constructed from ladder operators for the associated Laguerre polynomials. 
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As can be observed from (5.26) and (5.27), the appropriate ladder operators for the 



wavefunctions can be constructed from ladder operators for the Laguerre polynomials 
that shift the degree by one and the parameter by two. To obtain such ladder operators, 
we take combinations of the following ladder and shift operators, see e.g. [I], 

d 



x— + [a - x 
ox 



a+l , 
n-V 



X 



d 

x— + a 
ox 



d_ 

dx 



- 1 



L a n (x) = (n + l)L a n -\> 
L a n (x) = (n + a)L a r r 1 ' 

m( X ) = -li + \x). 



(5.28) 
(5.29) 

(5.30) 

(5.31) 



Taking (5.28) composed with (5.31) and ( |5.29 ) composed with (5.30), modulo the 
eigenvalue equation for the Laguerre polynomials 



' d 2 M , d 

+ {^- + a- x)— + n 
ox 1 ox 



L a Jx) 



0. 



(5.32) 



yields first-order operators with the appropriate action on the Laguerre polynomials. 
Conjugating by an appropriate factor, gives the following ladder operators 



, d 2n + a + 1 a(l + a) 

1 + 0)7— H 

'dx 2 2x 

d 2n + a + 1 a(a — 1) 

~ a '~dx~ + 2 2x 



e-^xH*(x) = -e-^x% +1 L a n t 2 1 (x) (5.33) 
e~$x%L*(x) = (n+l)(n+a)e~*x*~ 1 L%? 1 (x). (5.34) 



For the function i? m ,„(£), the ladder operators (5.33), (5.34) take the following 
forms 



K 



V-A n ,E„ 



V /Z A^(1 + V^Ar) (5.35) 



K 



V~An,E n 



= (l - - ^ + " V^)- (5-36) 



1 

2? 



The action of the K's on the functions R m n is as follows 



~ Rm—l,n— 1/k 

Thus, to raise or lower the index m by p, we need only to apply the operator p times 



K_.r^v R ill .1 



K p 



K p 



V—A n ,E n 



K 



y/^A^~2(p-l),E„ 



^ ' —^/~A n — 2,E m . n ~ i s/—%An,Er ni ,ti 



Note that after each successive application, we change the operator by adding 2 to 
\J — A n but the quantity E m>n is unchanged in each iteration. Note that by adding or 
subtracting p from a/— A n we change n by q 



y^A~ n ±2p = P - \-2n + - - 1 ) ± 2p = yJ-A nTq , 

q \ 2y/—ak J 



(5.37) 
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and so the action of these operators K p is given by 

^^/^A^,E m n -^ m .™ = ( _ ^-Y Rm+p,n+q (5.38) 

K P , — j — Rm.n = ( m 4" l)p(l m V A n ) pR m _„ n -q ■ (5.39) 

Similarly, we have raising and lowering operators for the functions S n (rj) given by 

Jv=K,k = (l + V^An/k) V - -JL^ - JL v^( fc + v^) (5-40) 

J-V^/k = (l - V^A) - + ^V 73 ^ - >/=40- (5-41) 
The action on the basis is 



Again, the repeated application of the ladder operators is defined as 



J -V^/k S n = ~{n + l)(-n + - l)5 n+ i. 



with action on as 

= (-l)^n-, (5.42) 

= ( n + ^ - 2^k )qS ^ (5 ' 43) 
Thus, we have created parameter dependent raising and lowering operators which 
map between wave fucntions with the same energy values given by 

To remove the dependence of these operators on the indices m, n, first we push the 
constant E m>n to the right and then replace with the operator H. Thus, when the 
resulting operator acts on a wave function ^ m ,n, it takes the value of the energy E m>n 
Next, we note that after this replacement, the resulting operators 

T s = Tx + T 2 , Ta = —L=(T l -T 2 ), (5.46) 

V~ A n 

are polynomials in A n (and hence the index n) and so we can again push this constant 
to the right and replace it with the first integral of the motion 

A n -> L x = -d 2 s - ae 4ks - (3e 2ks 

to remove the dependence of the operator on n. In this way, we have removed the 
dependence of the operators on the quantum numbers m and n and so 

[T s ,H k }^ min = [T a ,H k }^ m>n = 0, Vm,n. (5.47) 
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By a standard arguments as given in [22] , the identities hold in general (not only on the 
solutions of the Schrddinger equation) 

[T s ,H k ] = [T a ,H k }=0, (5.48) 

Thus, we have proven that the system is superintegrable. 

6. Conclusions 

In Section 2 we have proposed a method for expanding an integrable and exactly 
solvable system that separates in subgroup type coordinates into an infinite family 
of such systems. In the article, we have shown that in two-dimensional flat spaces the 
expanded family remains in the same space (E 2 (C), E 2 or £1,1) and that the method also 
preserves superintegrability. The method can also be used in higher-dimensional spaces 
and in non-flat ones. It will always embed a given integrable and exactly solvable system 
into a parameterized family of such systems. In general, it will however not preserve the 
structure of the space, i.e. it can change the metric. There is also no guarantee that when 
applied to a given superintegrable system, the method will preserve superintegrability. 

In Sections 3 and 4, the superintegrability of the classical and quantum system 
is demonstrated. For the classical system, the additional integral of the motion was 
constructed from the bounded trajectories thus proving the converse of the theorem of 
Nekhoroshev [36] for this system. For the quantum system, the exact-solvability of the 
wave- functions was used to construct the additional integral, thus proving the converse 
of the conjecture of Tempesta, Turbiner and Winternitz [44J for this system. 



The exact-solvability of the quantum system (5.1) in terms of generalized 
Bessel polynomials, provides yet another instance of the deep connection between 
superintegrable systems and orthogonal polynomials. Beyond the classical families of 
orthogonal polynomials known to be associated with superintegrable systems, there 
have been recent advances in discovering the connection between superintegrable 
systems and new families of orthogonal polynomials. In particular, new families of 
superintegrable systems associated with the q — > — 1 limit of q-Jacobi polynomials 
[4*0] and with exceptional Jacobi polynomials The connection between these 

two fields is evidenced in the methods of constructing the integrals of the motion 
from the ladder operators for the polynomials, as above and in [22l [33] 134"] . as well 
as in the representation theory of the algebras generated by such integrals see e.g. 
[18] [38] . Expanding the connection even further, are recent results linking other special 
functions, including elliptic and Painleve functions, to superintegrable systems, see e.g. 

[121 H31 EHl [321 EH H7]. 



Appendix A. Trajectories of the harmonic oscillator in Ei^ 



We give the trajectories of the harmonic oscillator in £7i,i, which were previously studied 
in [37]. In Cartesian coordinates the Hamiltonian is simply the difference of two one- 
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dimensional oscillators 

tjO 2,22222 / a i \ 

H =P u + oju — p v — uv. (A.l) 

The trajectories, up to rotation of the plane and translations in time, are given by 

u(t) = asm2out, v(t) = bcos2ut, (A. 2) 

p u (t) = uacos2ut, Pv(t) = —bu sin2ut. (A. 3) 

The energy can take positive and negative values and is defined by 

E = u 2 (a 2 -b 2 ). (A.4) 

In the coordinate system (p, a) the trajectories are given by 

p(t) = a 2 sin 2 2cot — b 2 cos 2u>t, 

, . a sin 2ut + b cos 2ut 
a(t) = 



P P (t) = 



a sin 2ut — b cos 2cut ' 
uo(a 2 + b 2 ) sin 2ut cos 2wt 
2 (a 2 sin 2 2ut - b 2 cos 2wt) 



, . uab(a sin 2wt — 6 cos 2ut) 
' ~ 2(asin2wt + 6cos2wt) ^ ' ^ 

with the separation constant A = 4(T 2 p 2 taking values A = u 2 a 2 b 2 . 
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